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ABSTRACT 



We present the first high-resolution X-ray spectrum of the Seyfert 1.5 galaxy NGC 
4151. Observations with the Chandra High Energy Transmission Grating Spectrometer 
reveal a spectrum dominated by narrow emission lines from a spatially resolved (1.6 
kpc), highly ionized nebula. The X-ray narrow- line region is composite, consisting of 
£f) . both photoionized and collisionally ionized components. The X-ray emission lines have 

similar velocities, widths, and spatial extent to the optical emission lines, showing that 



they arise in the same region. The clouds in the narrow-line region must contain a large 



range of ionization states in order to explain both the optical and X-ray photoionized 



emission. Chandra data give the first direct evidence of X-ray line emission from a hot 
plasma (T ~ 10 7 K) which may provide pressure confinement for the cooler (T = 3 x 10 4 
K) photoionized clouds. 

Subject headings: galaxies: individual (NGC 4151) — galaxies: Seyfert — X-rays: galaxies 



2 1 K) photoionized clouds. 

CO 



1. Introduction 

NGC 4151 (z=0. 00332) is a well-studied Seyfert 1.5 galaxy. It is known to be extended in soft 
X-rays from Einstein and ROSAT observations (Elvis, Briel, & Henry 1983; Morse et al. 1995). 
The nature of the extended X-ray emission was postulated to be either thermal emission from a 
collisionally ionized plasma or scattered nuclear flux. Weaver et al. (1994) modeled the ASCA 
spectrum by 3 components: an absorbed power law plus reprocessed emission from a warm absorber 
and an extra scattered component. Variability was observed in the column density of the warm 
absorber. 

Spectral complexity from unresolved emission lines and ionized absorption edges has been seen 
in the residuals of power-law fits to the spectra of many Seyfert galaxies. Recent observations with 
Chandra have revealed narrow absorption and emission lines in the spectra of Seyfert 1 galaxies 
NGC 5548 (Kaastra et al. 2000) and NGC 3783 (Kaspi et al. 2000). These two sources are 
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dominated by continuum emission from the central source, so absorption lines are more prominent 
than emission lines. On the other hand, the nucleus of NGC 4151 is highly absorbed, giving a good 
view of its narrow emission line region. 

In this paper we present Chandra High Energy Transmission Grating Spectrometer (HETGS) 
observations of NGC 4151. We identify and measure a large number of narrow emission lines in 
the X-ray spectrum. The high spectral resolution provided by HETGS allows us to study the 
ionization, temperature, and kinematics of the X-ray emitting regions. We discuss the association 
of the extended X-ray emission with the narrow-line region, and interpretations of the multiple 
X-ray emitting components in the spectrum. 

2. Observations 

We observed NGC 4151 for 48 ks on 2000 March 5 with Chandra HETGS (Canizares et al. 
2000), at a roll angle of 153.2°. The dispersion directions of the High Energy Grating (HEG) 
and the Medium Energy Grating (MEG) made acute angles of 42.5° and 37.5°, respectively, with 
the direction of the extended X-ray emission. Since the source has appreciable extent along the 
dispersion direction (see below), the spectral and spatial dimensions are mixed. Spectral extractions 
and reductions were performed using CIAO 1 and IDL. We used extraction windows of 7'.'9 width 
in the cross-dispersion direction, which were centered on the nucleus of the 0-order image. The 
windows included most of the flux from the extended source. We then summed the plus and minus 
orders of the spectra. The flux calibration of HETGS is currently accurate to 10% for energies 
E > 2 keV and 20% for E < 2 keV. 

The wavelength calibration was tied to an HETGS observation of Capella which was reduced in 
an identical manner to NGC 4151 (using CIAO, v. 2.0a, 2000 August 6). In this way we determined 
that velocities must be corrected by — 210±20 km s _1 to remove a systematic wavelength calibration 
error that exists in standard-processed data. This error is currently attributed to an error in the 
CCD pixel scale, due to thermal contraction 2 . Corrected velocities are reported in Table 1. 
Marginally redshifted X-ray narrow emission lines have been reported for NGC 5548 (Kaastra et 
al. 2000, 270 ± 100 km s" 1 ) and NGC 3783 (Kaspi et al. 2000, 230 ± 170 km s" 1 ). These redshifts 
may be owing to similar wavelength calibration errors. 

NGC 4151 was in a low flux state during our observation, with 3 x 10~ 3 photons s _1 cm~ 2 
A -1 at 2 A. We show the first order HEG and MEG spectra in Figure 1. Above 6 A (below 2.1 
keV), the spectrum is dominated by emission lines. Table 1 gives a list of the identified emission 
lines, their fluxes, and velocities relative to the galaxy rest frame. 



1 http://asc. harvard.edu/ciao 
2 http://asc. harvard.edu/cal 
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We compare the HETGS 0-order X-ray (0.4-2.5 keV) image of NGC 4151 to an optical [O 
in] 5007 image from the Hubble Space Telescope (HST) (Fig. 2). There is a good, but not exact 
correspondence between the structure of the optical and soft X-ray emission regions. The X-ray 
emission has a similar extent (23" = 1.6 kpc) and orientation (PA= 63°) to the optical narrow- 
line region (NLR). Most (70%) of the soft X-ray emission is resolved by Chandra; the rest comes 
from an FWHM= 0'.'9 nucleus. The hard X-ray emission (E>2.5 keV) is consistent with a nuclear 
point source. 

The spatial profiles of the narrow emission lines are extended (Figure 3). We show three lines, 
from three different ionization states: O vn (f), O vm Lya, and Fe I Ka. All three lines have 
similar spatial profiles to the 0-order soft X-ray image, and can be traced out to a distance of 
3" = 210 hjQ pc from the nucleus. 

HETGS is rather insensitive to relativistically broadened Fe Ka emission. Yaqoob et al. 
(1995) fit ASCA observations of Fe Ka with broad {a = 0.7 keV) and narrow components having 
equivalent widths of ~ 300 eV and ~ 100 eV, respectively. We find an equivalent width of 160 ± 20 
eV for narrow Fe I Ka with the HEG. The narrow core of the line is unresolved, with FWHM= 
1800 ± 200 km/s. 

3. Discussion 

3.1. Photoionized X-ray Line Emission 

The spectrum of NGC 4151 is hybrid, with line emission coming from both photoionized and 
collisionally ionized plasmas. There are several strong indications of a photoionized component. 
First, we see narrow radiative recombination continua (RRCs) from N vn, O vn, O vm, Ne ix, 
and Ne X. We estimate half widths (at 37% peak flux) for the RRCs which range from 1-4 eV, 
corresponding to temperatures of 2 — 4 x 10 4 K. To have such a low temperature and contain such 
high ionization states, the plasma where the RRC originates must be dominated by photoionization. 

The ratio of the O vn forbidden to resonance emission lines (f/r = 4.6 ± 1.5) is also consistent 
with photoionized plasma (Porquet & Dubau 2000). In addition, relatively strong n = 3 — 1 
transitions of N VI, Ne IX, and Mg XI are indicative of photoionization (Bautista & Kallman 
2000). The strength of O, Ne, Mg, and Si K-shell emission lines and relative weakness of Fe L-shell 
lines are additional characteristics of a photoionized plasma (Liedahl 1999). 

3.2. Thermal X-ray Line Emission 

Except for O vn, the He-like lines have a small f/r ratio, which means that collisional ionization 
is important. The Ne ix lines have f/r = 1.7±0.4. A ratio f/r ~ 3.4 is expected for a photoionized 
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plasma, and f/r < 0.8 for a collisionally ionized plasma with T > 2 x 10 6 K. There is clearly a 
mixture of line emission from both types of plasma. We note that fluorescent excitation by the 
continuum and resonance scattering may lead to enhancement of resonance lines under special 
circumstances (Bautista et al. 1998; Kallman et al. 1996). However, we would expect these 
processes to equally effect the O vn and other He-like resonance lines, which is not the case. 

The Lya lines of O VIII, Ne x, and Mg xn also have a large contribution from collisionally 
ionized plasma. The O vm Lya line is much stronger than the O vm RRC, while a ratio near 
unity is expected for a pure photoionized plasma. We find that 89 ± 5% of the O vm Lya emission 
comes from the hot plasma component of the NLR. Similarly, 70 ± 10% of the Ne x Lya emission 
comes from this component. A large fraction of the Mg xn Lya flux must also have a collisional 
origin, since the Mg xn RRC (6.32 A) is not even detected. 

From the range of formation temperatures of the collisionally augmented lines, we deduce 
that there is plasma in the temperature range T = 0.3 - 1 x 10 7 K. The nondetection of S XVI 
emission (4.73 A, 2 ± 5 x 10 -6 photons cm" 2 s" 1 ) is partly owing to the strong continuum at that 
wavelength. However, the detection of Fe xxv (Tf = 8 x 10 7 K) suggests that hotter gas is present. 
A more thorough analysis will be necessary to derive the volume emission measure vs. temperature 
distribution of the hot plasma. 

3.3. Continuum Emission 

We characterize the continuum emission from NGC 4151 with the sum of two absorbed power 
laws. The hard continuum (E > 3 keV) has a photon index T = 0.4 ± 0.3 and normalization 
A = 2.3^g'g x 10" 3 photons keV" 1 cm" 2 s" 1 at 1 keV, absorbed by a neutral column of Njj = 
3.7^o g x 10 22 cm" 2 . The observed 2-10 keV flux is 5.5 x 10" 11 ergs cm" 2 ,?" 1 . The soft continuum 
(E < 1.3 keV) has a photon index of T = 3.1 ± 0.5 and normalization A = 1.0 ± 0.1 x 10" 3 photons 
keV -1 cm" 2 s" 1 at 1 keV, absorbed by a Galactic column of Nh = 2.0 x 10 20 cm" 2 . 

We deduce that there is strong absorption of the hard nuclear flux by a cloud which does not 
occult the bulk of the soft X-ray emitting NLR. Weaver et al. (1994) observed a 2-10 keV flux 
of 11 - 22 x 10" 11 ergs cm"^" 1 , T ~ 1.5, and N H = 2 - 5 x 10 22 cm" 2 with ASCA. The hard 
continuum was harder and a factor of 2 — 4 times fainter during our observation, but Njj was 
similar. However, we note that poor statistics for energies above 7 keV in our data may bias the 
fit value for T. Weaver et al. (1994) also found variations in Nh with ASCA on a time scale of 
6 months, which suggests that absorption is taking place on a spatial scale much smaller than the 
NLR. 

The soft continuum of NGC 4151 has a much steeper photon index than the hard continuum 
(3.1 vs 0.4). This argues against Thomson scattered emission as the source of the soft continuum, 
unless the intrinsic nuclear continuum steepens drastically below 3 keV. Thermal bremsstrahlung 
and blended Fe L emission lines are therefore the most likely sources of the soft continuum. 
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3.4. Kinematics of the X-ray Line Region 

The locations and velocity field of optically emitting clouds in the NLR of NGC 4151 have 
been mapped in detail by Kaiser et al. (2000), and are consistent with a biconical outflow from 
the nucleus. Integrated long-slit spectroscopy of the central 2" of the NLR gives a blueshift of -90 
km s" 1 and a line width (FWHM) of 460 km s" 1 for the [O in] 5007 line (Heckman et al. 1981). 
The blueshifts and blue asymmetry of the optical lines are likely owing to obscuration of the far 
cone by dust. 

The unblended X-ray emission lines have an (unweighted) mean blueshift of —120 ±50 km s _1 
with respect to the galaxy rest frame, as measured from its stellar spectrum. This is consistent 
with the velocity of the optical NLR. Since the X-ray line widths are affected by spatial broadening, 
we can only get an upper limit to their velocity widths from the widths of the low energy lines. O 
vii (f) has an apparent width (Gaussian FWHM) of 460 km s _1 , which equals the widths of the 
optical lines. 

It is interesting that the X-ray emission region is more extended in the front (SW) cone than 
the back cone (Fig. 2), suggesting that the soft X-rays from the back cone are obscured, like the 
optical emission. The similar spatial extent of the X-ray and optical emission line regions (Figs. 
2,3) and their similar kinematics are strong evidence that they are closely linked. X-ray emission 
from the NLR demonstrates that high energy processes are important in this region, at distances 
of 70-800 pc from the galactic nucleus. 

3.5. The Multiphase Nature of the NLR 

For the first time, there is direct evidence for at least 2 distinct gas phases in the extended 
NLR of a Seyfert galaxy. One of the phases is cool (T = 3 x 10 4 K) and photoionized, while the 
second is hot (T ~ 10 7 K) and collisionally ionized. In addition, there is a large range in ionization 
of the photoionized region, as demonstrated by O i-O in in the optical-UV spectrum and O VII in 
the X-ray spectrum. There are also strong Ka lines from neutral species in the X-ray spectrum, 
including Fe I, Si I, and Mg I. Detailed modeling will be necessary to determine the volume emission 
measure vs. ionization distribution of the photoionized component. A wide range of ionization is 
derived from photoionization models of the NLR of the Circinus galaxy (Sako et al. 2000), and 
might be explained by a radial gradient in the NLR density. 

A large percentage (65 ± 9%) of the narrow Fe I Ka emission comes from the extended NLR 
(Fig. 3). This is contrary to the popular idea that the narrow emission comes primarily from a 
parsec scale torus. The remaining unresolved emission, which has an equivalent width of 56 ± 9 eV, 
may very well come from a torus. Unlike NGC 5548 (Yaqoob et al. 2000), there is no evidence for 
Doppler broadening of the narrow Fe I Ka line in NGC 4151. 

It has been suggested that a hot, inter-cloud medium is necessary for pressure confinement of 
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the NLR clouds (Elvis, Briel, & Henry 1983). The hot plasma we detect in our X-ray spectrum 
may form this medium. We derive volume emission measures for the thermal contribution to the 
Ne x and O vm Lya lines of n 2 e V = 2 — 4 x 10 63 /iy 2 cm -3 , assuming T = 1 x 10 7 K, and line 
emissivities for solar abundances from Mewe, Gronenschild, & van den Oord (1985). Crenshaw 
et al. (2000) derive a half-opening angle of 36° from kinematic models of the optical NLR. If the 
thermal X-ray line emission comes from a uniformly filled (/ = 1) bicone in the central 4" of the 
galaxy, then the mean plasma density is n e = 3 cm -3 and the pressure is n e T = 3 x 10 7 cm -3 
K. Pressure equilibrium between the hot and cool phases of the X-ray NLR would then require a 
density of n e = 10 3 cm' 3 for the cool phase. 

The optical NLR clouds at r = 6 — 20" have a density n e = 220 cm -3 and a pressure of 
n e T = 3 x 10 6 cm -3 K (Penston et al. 1990; Morse et al. 1995), much lower than what we 
find closer to the nucleus. However, we know that the spatial profile of the NLR flux is centrally 
peaked. A density profile of r -2 would account for the larger pressure of the nucleus than the 
extended NLR, assuming constant temperature. Estimates of the pressure can be improved by 
modeling the density distribution to match the spatial profiles of the X-ray emission. 

A possible alternative is that the hot plasma is confined to shocks between the NL clouds and 
the host galaxy ISM. In that case, its filling factor / would be much lower, and its density and 
pressure much higher than calculated above. The X-ray emission we observe in the Chandra image 
has a clumpy morphology which follows the optically emitting clouds. This may be partly owing 
to emission from a shocked component. There is evidence for deceleration of the NL clouds at a 
radius of about 2" from the nucleus (Crenshaw et al. 2000), perhaps from interaction with the 
ISM. 

4. Conclusions 

The power of Chandra HETGS to elucidate the nature of the X-ray emission from active 
galactic nuclei is demonstrated by high resolution spectra and images of NGC 4151. Direct emission 
from the nucleus is highly absorbed. Strong, narrow X-ray emission lines are seen from the spatially 
resolved NLR of NGC 4151. The narrow radiative recombination continuum features and the O 
vii f/r ratio are consistent with photoionization from the hard nuclear X-ray source. There is 
also narrow line emission from a hot, collisionally ionized component in the extended NLR. This 
may come from the long-sought intercloud medium or shocked plasma. In addition, we find strong 
narrow Fe I Ka emission from the extended NLR. The composite nature of the NLR indicates that 
both photoionization and collisional heating are important. The extent and kinematics of both 
X-ray photoionized clouds and hot thermal plasma match those of the optical NLR, showing that 
they are closely connected. We also suggest that the narrow X-ray emission lines seen in some 
other Seyfert galaxies may come from the NLR, and may be interpreted in the light of these NGC 
4151 observations. 
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Table 1. NGC 4151 Emission Lines 



A(rest A) 


v(km s 1 ) 


Flux a 


SNR 


b ID 


1.79 


460±420 


3.3E-05 


2.3 


Fe xxv 


1.937 


250±80 


1.8E-04 


11.0 


Fe I Ka 


6.182 


-20±120 


1.0E-05 


4.1 


Si xiv Lya 


6.648 


-30±130 


1.1E-05 


4.4 


Si xni r 


6.740 


-40±40 


1.1E-05 


5.0 


Si xni f 


7.106 


-380±110 


5.3E-06 


3.5 


Mg xn Ly/3 c 


7.130 


-380±40 


1.3E-05 


5.6 


Si I Ka 


7.851 


-560±80 


2.5E-06 


2.4 


Mg XI ls3p-ls 2 


8.421 


-210±90 


1.1E-05 


7.4 


Mg xn Lya 


9.01 


-380±170 


3.5E-06 


3.3 


Fe xxn c 


9.102 




7.2E-06 


3.2 


Ne x RRC c 


9.169 


-20±160 


9.3E-06 


6.3 


Mg xi r c 


9.314 


-406±140 


1.2E-05 


7.4 


Mg xi f 


9.887 


-330±50 


3.2E-06 


2.7 


Mg I Ka 


10.239 


-180±70 


7.4E-06 


4.5 


Ne x Ly/3 


10.368 




9.3E-06 


4.2 


Ne ix RRC 


11.001 


110±140 


3.8E-06 


2.2 


Ne ix ls4p-ls 2 


11.547 


-290±70 


5.9E-06 


3.0 


Ne ix ls3p-ls 2 


12.134 


-70±80 


2.0E-05 


6.4 


Ne x Lya 


13.447 


-150±40 


2.0E-05 


4.8 


Ne ix r 


13.553 


-210±160 


9.9E-06 


2.5 


Ne ix i 


13.698 


10±40 


3.3E-05 


6.1 


Ne ix f 


14.228 




1.1E-05 


2.3 


vin RRC c 


15.176 


-170±50 


1.1E-05 


3.2 


viii Ly7 


16.006 


-250±20 


1.8E-05 


3.8 


viii Ly/3 c 


16.771 




5.2E-05 


4.7 


vii RRC 


18.588 




2.6E-05 


2.6 


N vii RRC 


18.969 


-40±30 


1.0E-04 


7.0 


viii Lya 


21.602 


106±20 


6.8E-05 


4.1 


vn r 


21.804 


-410±140 


5.1E-05 


3.3 


vii i 


22.101 


-132±20 


3.1E-04 


8.3 


vii f 


24.781 


-190±60 


6.3E-05 


4.0 


N vii Lya 


24.898 


-10±20 


3.9E-05 


3.0 


N vi ls3p-ls 2 



a Flux (photons s 1 cm 2 ) 

b Signal-to- noise ratio 

c Line blend or possible line blend 
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Figure Captions 

Figure 1 

Chandra first order HEG (top) and MEG (middle, bottom) spectra of NGC 4151, with line identi- 
fications. Above 6 A, the spectrum is dominated by narrow emission lines from both thermal and 
photoionized plasmas. Note that the ratio of forbidden to resonance emission in the He- like O vn 
lines is large, while the ratio for Ne ix, Mg XI, and Si xiii is smaller, indicating a hybrid spectrum. 
The narrow radiative recombination continua (RRCs) of several ions attest to photoionization in a 
cool T = 3 x W^K plasma. 

Figure 2 

Comparison of X-ray and optical emission from the extended NLR in NGC 4151. Chandra contours 
(0.4-2.5 keV band) are overlayed on an HST [O m] 5007 image. A 36'/ 4 by 23'.'5 region is shown, 
with N at top and E at left. Note that the SW (front) cone is more extended than the NE (back) 
cone in both bands. X-ray contours are separated by factors of 2 in surface brightness, with the 
lowest contour at 2.1 photons arcsec -2 . 

Figure 3 

Comparison of NGC 4151 zero-order spatial profile and 1st order emission line spatial profiles. The 
gaussian profile of the unresolved nuclear component is shown for comparison. The spatial profiles 
of O vii, O viii, and Fe Ka are all resolved, from a region at least 6" in extent. 
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Fig. 1. — See figure caption. 
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Fig. 3. — See figure caption. 



